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Abstract. A few characteristic exponents describing power law behaviors of rough- 
ness, coherence and persistence in stochastic time series are compared to each other. 
Relevant techniques for analyzing such time series are recalled in order to distin- 
guish how the various exponents are measured, and what basic differences exist 
between each one. Financial time series, like the JPY/DEM and USD/DEM ex- 
change rates are used for illustration, but mathematical ones, like (fractional or 
not) Brownian walks can be used also as indicated. 



1 Introduction 

A great challenge in modern times is the construction of predictive theories 
for nonlinear dynamical systems for which the evolution equations are barely 
known, if known at all. General or so-called universal laws are aimed at from 
very noisy data. A universal law should hold for different systems charac- 
terized by different models, but leading to similar basic parameters, like the 
critical exponents, depending only on the dimensionality of the system and 
the number of components of the order parameter. This in fine leads to a 
predictive value or power of the universal laws. 

In order to obtain universal laws in stochastic systems one has to distin- 
guish true noise from chaotic behavior, and sort out coherent sequences from 
random ones in experimentally obtained signals [Q. The stochastic aspects 
are not only found in the statistical distribution of underlying frequencies 
characterizing the Fourier transform of the signal, but also in the amplitude 
fluctuation distribution and high moments or correlation functions. 

Following the scaling hypothesis idea, neither time nor length scales have 
to be considered Henceforth the fractal geometry is a perfect frame- 

work for studies of stochastic systems which do not appear at first to have 
underlying scales. A universal law can be a so-called scaling law if a log[y{x)] 
vs. log{x) plot gives a straight line (over several decades if possible) leading 
to a slope measurement and the exponent characterizing the power law. 

When examining such phenomena, it is often recognized that some coher- 
ent factor is implied. Yet there are states which cannot be reached without 
going through intricate evolutions, implying concepts like transience and per- 
sistence, - well known if one recalls the turbulence phenomenon and its basic 
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theoretical understanding g. Finally, the apparent roughness of the signal 
can be put in mathematical terms. These concepts are briefly elaborated 
upon in Sect. 2. 

In the y{x) function, y and x can be many "things" . However the relevant 
outlined concepts can be well illustrated when x is the time t variable. In so 
doing time series serve as fundamental testing grounds. Several series can be 
found in the literature. Financial time series and mathematical ones based on 
the Weierstrass-Mandelbrot function Q describing fractional (or not) Brow- 
nian walks can be used for illustration. The number of points should be large 
enough to obtain small error bars. A few useful references, among many oth- 
ers, discussing tests and other basic or technical considerations on non linear 
time series analysis are to be found in @,^,^,^o|jl^,|l^ . 

Mathematical series, like (fractional or not) Brownian motions (Bm) and 
practical ones, like financial time series, are thus of interest for discussion or 
illustrations as done in Sect. 3. There are several papers and books for interest 
geared at financial time series analysis ... and forecasting. Again not all can 
be mentioned, though see ||l^,^^Jl^. On a more general basis, an introduction 
to financial market analysis per se, can be found in JlGf . 

Nevertheless it should be considered that any "scaling exponent" should 
be robust in a statistical sense with respect to small changes in the data or 
in the data analysis technique. If this is so some physics considerations and 
modeling can be pursued. One question is often raised for statistical pur- 
poses whether the data is stationary or not, i.e. whether the analyzed raw 
signal, or any of its combinations depend on the (time) origin of the series. 
This theoretical question seems somewhat practically irrelevant in financial, 
meteorological, ... sciences because the data is obviously never stationary. In 
fact, in such new "exotic applications of physics" a restricted criterion for 
stationarity is thought to be sufficient: if the data statistical mean and the 
whatever-extracted-parameter do not change too much (up to some statis- 
tical significance Q) the data is called quasi — stationary. If so it can be 
next offered for fundamental investigations. Thereafter, the prefix "quasi" is 
immediately forgotten and not written anymore. 

Several characteristics plots leading to universality considerations, thus 
fractal-like exponents are first to be recalled. They are obtained from different 
techniques which are briefly reviewed for completeness either in Sect. 3 or in 
Appendices. Some technical materials can be usefully found in This 
should serve to distinguish how the various exponents are measured, and what 
basic differences exist between each one. The numerical values pertaining to 
the words (i) persistence, (ii) coherence, and (iii) roughness will be given and 
related to each other. From a practical point of view, as illustrated in the 
exercice session which was taking place after the lecture, the cases of foreign 
exchange currency rates, i.e. DEM/USD and DEM/JPY are used. They are 
shown in Fig.l and 2 for a time interval ranging from Jan. 01, 1993 till June 
30, 2000. 
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Fig. 1. Evolution of the exchange rate DEM/USD from Jan. 01, 1993 till June 30, 
2000 
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Fig. 2. Evolution of the exchange rate DEM/JPY from Jan. 01, 1993 till June 30, 
2000 



For an adequate perspective, let it be recalled that such modern con- 
cepts of statistical physics have been recently applied in analyzing time series 
outside finance as well, like in particular those arising from biology |2C,pT|, 
medicine |^^, meteorology |2^, electronics |Q, image recognition |25|, ... 
again without intending to list all references of interest as should be done 
in a (longer) review paper. Many examples can also be found in this book 
through contributions by world specialists of computer simulations. 
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2 Phase, amplitude and frequency revisited 
2.1 Coherence 

When mentioning the word coherence to any student in physics, he/she is 
immediately thinking about " lasers" , - that is where the word has been most 
striking in any scientific memory. It is recalled that the laser is a so interest- 
ing instrument because all photons are emitted in phase, more exactly the 
difference in phases between emitted waves is a constant in time. Thus there 
is a so-called coherence in the light beam. They are other cases in which 
phase coherence occurs, let it be recalled that light bugs are emitting coher- 
ently [ p6| , young girls in dormitories have their period in a coherent way, 
driving conditions are best if some coherence is imposed |^ , crystals have a 
better shape and properties if they are grown in a " coherent way" ; sand piles 
and stock markets seem also to have coherent properties pq,E9|. 



2.2 Roughness 

No need to say that a wave is characterized by its amplitude which has also 
some known importance in measurements indeed, be it often the measure 
of an intensity (the square of the amplitude). More interestingly one can 
define the roughness of a profile by observing how the signal amplitude varies 
in time (and space if necessary), in particular the correlation between the 
various amplitude fluctuations. 



2.3 Persistence 

On the other hand, it can be easily shown that a periodic signal can be de- 
composed into a series of sin or cos for which the frequencies are in arithmetic 
order. Thus, the third "parameter" of the wave is its frequency. Some (reg- 
ular) frequency effect is surely apparent in all cycling phenomena, starting 
from biology, climatology, meteorology, astronomy, but also stock markets, 
foreign currency exchange markets, tectonics events, traffic and turbulence, 
and politics. For non periodic signals, the Fourier transform has been intro- 
duced in order to sort out the distribution of frequencies of interest, i.e. the 
" density of modes" . The distribution defines the sort of persistence of a phe- 
nomenon. It might be also examined whether the frequencies are distributed 
in a geometrical progression, rather than following an ordinary/usual arith- 
metic progression, i.e. whether the phenomena might be log-periodical, like 
in antennas, earthquakes and stock market crashes ]30|,^l||. 



3 Power Law Exponents 
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3.1 Persistence and Spectral Density 

Data from a (usually discrete) time series y(t) are one dimensional sets and 



are more simple to analyze at first than spatial ones 33|. Below and for 
simplicity one considers that the measurements are taken at equal time inter- 
vals. Thus for financial time series, there is no holiday nor week-end. A more 
general situation is hardly necessary here. Two classic examples of a math- 
ematical univariate stochastic time series are those resulting from Brownian 
motion and Levy walk cases [Q. In both cases, the power spectral density 
Si{f) of the (supposed to be self-affine) time series y{t) has a single power-law 
dependence on the frequency /, 

^i(/) (1) 

following from the Fourier transform 

Si{f)= J dte^f'yit). (2) 
For y{t) one can use the Weierstrass-Mandelbrot (fractal) function Q 

W{t) = 7^"^^^"' [1 - e'^"*] e"^''*, (3) 

m 

with 7 > 1, and 1 < D < 2. The phase 0„ can be stochastic or determin- 
istic. For illustration, Berry and Lewis |^ took (f)n — n fi, with fi — Q or tt. 
The function obeys 

W(7i) = 7(2-^) e-'^" W{t), (4) 
is stationary, and its trend, in the deterministic cases, is given by 

Wf{t) - W/f{t) - t^-^/ln 7. (5) 
The power spectrum is easily calculated Q to be 

S,{f) ^ ...{l/ln 7) f"-". (6) 



One could search whether the function moments obey power laws with 
characterizing exponents. Equation(l) allows one to put the phenomena into 
the self-affine class of persistent phenomena characterized by the [3 value. 
The range over which (3 is well defined in Eq.(l) indicates the range of the 
persistence in the time series. A Brownian motion is characterized by /? = 2, 
and a white noise by ^ = 0.|3 See a very interesting set of such mathematical 
signals and the corresponding power spectrum in p2[ . 

^ Notice that the differences between adjacent values of a Brownian motion ampli- 
tude result in white noise. 



6 Marcel Ausloos 



DEM/USD 

10^ , ^ 




10-" 10-' 10-' 10-' 10° 

f(day-^) 

Fig. 3. The power spectrum of the DEM/USD exchange rate for the time interval 
data in Fig.l 
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Fig. 4. The power spectrum of the DEM/JPY exchange rate for the time interval 
data in Fig.2 
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The Fourier transform, or power spectrum, of the financial signals used for 
illustration here are found in Fig. 3 and Fig. 4. Notice the large error bars, 
allowing to estimate that (3 is about equal to 2, as for a trivial Brownian 
motion case. However the coherence and/or roughness aspect are masked in 
this one-shot analysis. Only the persistence behavior is touched upon. 

If the distribution of fluctuations is not a power law, or if marked devi- 
ations exist, say the statistical correlation coefficient is less than 0.99 , indi- 
cating that a mere power law for S{f) is doubtful, a more thorough search of 
the basic frequencies is in order. A crucial step is to extract deterministic or 
stochastic components Q, e.g. the stochastic aspects found in the statistical 
distribution of values show its persistence to be either nonexistent (white 
noise case) or existent, i.e. /? 2. If so the persistence can be qualitatively 
thought to be strong or weak. The Fourier transform can sometimes indi- 
cate the presence of specific frequencies, much more abundant than others, 
in particular if cycles exist, as in meteorology and climatology. 

The range of the persistence is obtained from the correlations between 
events. A "short" or "long range" is checked through the autocorrelation 
function, usually ci. This function is a particular case of the so-called "q—th 
order structure function" |]35[] or "g— th order height-height correlation func- 
tion" of the (normalized) time-dependent signal y(ti), 



where only non-zero terms are considered in the average (.)^ taken over all 
couples (ti+r, ti) such that r = \ti^r —ti\- In so doing one can obtain a set of 
exponents /3g 

If the autocorrelation is larger than unity for some long time t one can 
talk about strong persistence, otherwise it is weak. This criterion defines the 
scale of time t in y{t) for which there is long or short (time) range persistence. 
Notice that the lower limit of the time scale is due to the discretization step, 
and this sets the highest frequency to be the inverse of twice the discretizatin 
interval. The upper limit is obvious. 

3.2 Roughness, Fractal Dimension, Hurst exponent 
and Detrended Fluctuation Analysis 

The fractal dimensionciterefAl,addison, falconer, roughness D is often used 
to characterize the roughness of profiles Several methods are used for 
measuring D, like the box counting method, not quite but are not quite 
efficient; many others are found in the literature as seen in and 
here below. For topologically one dimensional systems, the fractal dimension 
D is related to the exponent (3 by 



c,{r)^Mt,+r)-y{t.)nj{\y{ti}n, 



(7) 



(3 ^5 -2D. 



(8) 
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A Brownian motion is characterized hy D — 3/2, and a white noise by D 
= 2.0|jl2|. 

Another " measure" of a signal roughness is sometimes given by the Hurst 
Hu exponent, first defined in the "rescale range theory" (of Hurst |3^,Q ) 
who measured the Nile flooding and drought amplitudes. The Hurst method 
consists in listing the differences between the observed value at a discrete 
time t over an interval with size N on which the mean has been taken. The 
upper {um) and lower (ym) values in that interval define the range Rn = 
Um — Um- The root mean square deviation Sn being also calculated, the 
"rescaled range" is Rn/Sn is expected to behave like iV^". This means that 
for a (discrete) self-affine signal y{t), the neighborhood of a particular point 
on the signal can be rescaled by a factor b using the roughness (or Hurst 
[HQ]) exponent Hu and defining the new signal b^^"y{bt). For the exponent 
value Hu, the frequency dependence of the signal so obtained should be 
undistinguishablc from the original one, i.e. y(t). 

The roughness (Hurst) exponent Hu can be calculated from the height- 
height correlation function ci (t) supposed to behave like 

ci(T) = (|y(n,+,)-y(n,)|), ~r^^ (9) 

whereas 

Hu=l+Hi, (10) 

rather than from the box counting method. For a persistent signal. Hi > 1/2; 
for an anti-persistent signal. Hi < 1/2. Flandrin has theoretically proved 
[!§ that 

(3^2Hu-l, (11) 

thus /? = 1+2 Hi. This implies that the classical random walk (Brownian 
motion) is such that Hu = 3/2. It is clear that 



D^i-Hu. (12) 

Fractional Brownian motion values are practically found to lie between 
1 and 2 |]l7| , p"8|j2^ . Since a white noise is a truly random process, it can be 
concluded that Hu = 1.5 implies an uncorrelated time series p^ . 

Thus D > 1.5, or Hu < 1.5 implies antipersistence and D < 1.5, or 
Hu > 1.5 implies persistence. From preimposed Hu values of a fractional 
Brownian motion series, it is found that the equality here above usually holds 
true in a very limited range and (3 only slowly converges toward the value 
Hu [||. 

The inertia axes of the 2-variability diagram seem to be related 

to these values and could be used for fast measurements as well. 

The above results can be compared to those obtained from the Detrended 
Fluctuation Analysis |||,||| (DFA) method. DFA ^ consists in dividing a 
random variable sequence y{n) over N points into N/t boxes, each containing 
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r points. The best Unear trend z(n) = an + b in each box is defined. The 
fluctuation function F(t) is then calculated following 

kr 

= - E \yin)-z{n)\\ k = 1,2, ■■■,N/t. (13) 

n=(fc-l)r+l 

Averaging F[t)^ over the N Jt intervals gives the fluctuations {F{t)'^) as a 
function of r. If the y{n) data are random uncorrelated variables or short 
range correlated variables, the behavior is expected to be a power law 

(F2)1/2 (14) 

with Ha different from 0.5. 

The exponent Ha is so-labelled for Hausdorff |^|4|. It is expected, not 
always proved as emphasized by that 

Ha = 2-D, (15) 
where D is the self-affine fractal dimension [p|j3^ . It is immediately seen that 

(3^1 + 2Ha. (16) 

For Brownian motion. Ha — 0.5, while for white noise Ha = and D — 2. 




10° io' io" lo" 10* 

t(days) 

Fig. 5. The DFA result for the DEM/USD exchange rate for the time interval data 
in Fig.l 

The DFA log-log plots of the DEM/JPY and DEM/USD exchange rates 
are given in Fig. 5 and Fig. 6. It is seen that the value of Ha fulfills the above 
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Fig. 6. The DFA result for the DEM/JPY exchange rate for the time interval data 
in Fig.2 



relations for the DEM/JPY and DEM/USD data, since Ha = 0.55. Both Ha 
and D readily measure the roughness and persistence strength. Fractional 
Brownian motion Ha values are found to lie between and 1. The effect of 
a trend is supposedly eliminated here. However only the linear or cubic de- 
trending have been studied to my knowledge [4^ . Other trends, emphasizing 
some characteristic frequency, like seasonal cycles, could be further studied. 
A generalized Hurst exponent H{q) is defined through the relation 

Cq(T) ocr9-f^(9), g>0 (17) 

where Cq{T) has been defined here above 

The intermittency of the signal can be studied through the so-called sin- 
gular measure analysis of the small-scale gradient field obtained from the 
data through 

e{r; I) = '-"^ ^^=r' ly(^»+r) - y(^»)l ^^g^ 

< \y(.ti+r) - y{ti)\ > 

i = 0,...,A-r (19) 

r = l,2,...,vl = 2™ , (20) 

where m is an integer. The scaling properties of the generating function are 
then searched for through the equation 

X,iT)^<s{r;iy>^r-"^'^\ q>0, (21) 



with 
and 
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with T as defined above. 

The iir(g)-exponent is closely related to the generalized dimensions Dq = 
1 — K{q)/{q — 1) The nonlinearity of both characteristics exponents, 
qH{q) and K{q), describes the multifractality of the signal. If a linear depen- 
dence is obtained, then the signal is monofractal or in other words, the data 
follows a simple scaling law for these values of q. Thus the exponent ||35|,p6| 



dq 



(22) 



is a measure of the intermittency lying in the signal y{n) and can be nu- 
merically estimated by measuring Kq around q = 1. Some conjecture on the 



role/meaning of Hi is found in |28 . From some financial and political data 
analysis it seems that Hi is a measure of the information entropy of the 
system. 



4 Conclusion 

It has been emphasized that to analyze stochastic time series, like those 
describing fractional Brownian motion and foreign currency exchange rates 
reduces to examining the distribution of and correlations between ampli- 
tudes, frequencies an phases of harmonic-like components of the signal. Due 
to some scaling hypothesis, characteristic exponents can be obtained to de- 
scribe power laws. The usefulness of such exponents serves in determining 
universality classes, and in fine building physical or algorithmic models. In 
the case of financial times series, the exponents can even serve into imagin- 
ing some investment strategy . Notice that due to the non stationarity of 
the data, such exponents vary with time, and multifractal concepts must be 
brought in at a refining stage, - including in an investment strategy. 

In that spirit, let it be emphasized here the analogy between a Ci di- 
agram and the w, k diagram of dynamical second order phase transitions |^ . 
In the latter the frequency and the phase of a time signal are considered on 
the same footing, and encompass the critical and hydrodynamical regions. In 
the present cases an analog diagram relates the roughness and intermittency. 
This has been already examined in p8| , ^ . 

Among other various physical data analysis techniques which have been 
recently presented in order to obtain some information on the deterministic 
and/or chaotic content of univariate data, let us point out the wavelet tech- 
nique which has been considerably used (several references exist, see below) 
including for DNA and meteorology studi es p6| , ^ol[5ll , |5^ . The Hi , Ci tech- 
nique used in turbulence and meteorology [p4| is also somewhat appropriate. 

There are many other techniques which are not mentioned here, like the 



weighted fixed point |55 , and the time-delay embedding M. Surely several 



others have been used, but only those relevant for the present purpose have 
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Table 1. Values of the most relevant exponents in various regimes (i.e., stationary, 
persistent, antipcrsistcnt) of univariate stochastic scries : D : fractal dimension; Ha 
: Hausdorff meEisure; Hu : Hurst exponent ; a : from DFA technique; /3 : power 
spectrum exponent: WN = white noise, {f)Bm : (fractional) Brownian motion; 
flat : flat spectrum. 



Signal name 


D 


Ha 


Hu 


a 


P 












0.5 


antipersist 
0.5 


-1 



station 
uncorrel 


WN 


2 





1 


persist 


1 


station 


fBm 
Bm 
flat 


3/2 
1 


0.5 
1 


3/2 
2 


superpersist 
superpersist 


2 
3 


nonstat 



been fully mentioned hereabove. As a summary of the above, Table 1 indi- 
cates the range of values found for different signals and their relationship to 
stationarity, persistence and coherence. In the many years to come, it seems 
relevant to ask for more data on multivariate functions, thus extending the 
above considerations to other real mathematical and physical cases in higher 
dimensions. 

Finally, two short Appendices should follow in order to remain consis- 
tent with the oral lectures. It was shown that another time series "analysis" 
technique is often used by experts for some predictability purpose, i.e. the 
moving average technique. It is briefly discussed in Appendix A. It has served 
in the iihungen. Also, these lecture notes would bo incomplete without men- 
tioning the intrinsic discrete scale invariance implication in time series. Such 
a substructure leads to log-periodic oscillations in the time series, whence 
to fascinating effects and surprises in predicting crash-like events. This is 
mentioned in Appendix B. 
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5 Appendix A : moving averages 

To take into account the trend can be shown to be irrelevant for such short 
range correlation events. The trend is anyway quite ill defined since it is a 
statistical mean, and thus depends on the size of the interval, i.e. the num- 
ber of data points which is taken into account. Nevertheless many technical 
analyses rely on signal averages over various time intervals, like the moving 
average method p7[ | . These methods shold be examined form a physical point 
of view. An interesting observation has resulted from checking the density of 
intersections of such mean values over different time interval windows, which 
are continuously shifted. This corresponds to obtaining a spectrum of the 



so-called moving averages |57 , used by analysts in order to point to "gold" 
or "death" crosses in a market. The density p of crossing points between any 
two moving averages is obviously a measure of long-range power-law correla- 
tions in the signal. It has been found that p is a symmetric function of AT ^ 
i.e. the difference between the interval sizes on which the averages are taken, 
and it has a simple power law form |38[|59|] . This leads to a very fast and 
rather reliable measure of the fractal dimension of th signal. The method 
can be easily implemented for obtaining the time evolution of D, thus for 
elementary investment strategies. 



6 Appendix B : discrete scale invariance 



It has been proposed that an economic index y{t) follows a complex power 
law [lom, i.e. 



y{t) ^ A + B {tc~t)-"' [I + C cos{uj ln{{tc-t)/tc) + 0)] (23) 

for t < tc, where tc is the crash-time or rupture point, A, B, m, C, uj, 
(j) are parameters. This index evolution is a power law (to) divergence (for 
m > 0) on which log-periodic [uj] oscillations are taking place. The law for 
y{t) diverges at t—tc with an exponent to (for to > 0) while the period of 
the oscillations converges to the rupture point at t — tc- This law is similar 
to that of critical points at so-called second order phase transitions , but 
generalizes the scaleless situation for cases in which discrete scale invariance 
is presupposed |6^ . This relationship was already proposed in order to fit 
experimental measurements of sound wave rate emissions prior to the rupture 
of heterogeneous composite stressed up to failure The same type of 
complex power law behavior has been observed as a precursor of the Kobe 
earthquake in Japan Q. 

Fits using Eq.(23) were performed on the S'&P500 data |^,^ for the 
period preceding the 1987 October crash. The parameter values have not 
been found to be robust against small perturbations, like a change in the 
phase of the signal. It is known in fact that a nonlinear seven parameter fit 
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is highly unstable from a numerical point of view. Indeed, eliminating the 
contribution of the oscillations in Eq.(l), i.e. setting C=0, implies that the 
best fit leads to an exponent m = 0.7 quite larger than m = 0.33 for C 7^ 
[ |60| . Feigenbaum and Freund [ |6T| also reported various values of m ranging 
from 0.53 to 0.06 for various indexes and events (upsurges and crashes). 

Universality in this case means that the value of m should be the same 
for any crash and for any index. In so doing a single model should describe 
the phase transition, and the exponent would define the model and be the 
only parameter. A limiting case of a power law behavior is the logarithmic 
behavior, corresponding to m = 0, i.e. the divergence of the index y for t 
close to tc should be 



y{t) = A + B ln{{tc-t)/tc) [1 + C cos{ujln{{tc-t)/tc) + (f>)] (24) 

This logarithmic behavior is known in physics as characterizing the spe- 
cific heat ("four point correlation function") of the Ising model, and the 
Kosterlitz-Thouless phase transition in spatial dimensions equal to two. 
They are thus specific to systems with a low order dimension of the order pa- 
rameter. It is nevertheless a smooth transition. The mean value of the order 
parameter is not defined over long range scales, but a phase transition 
nevertheless exists because there is some ordered state on small scales. In 
addition to the physical interpretation of the latter relationship, the advan- 
tages are that (i) the number of parameters is reduced by one, and (ii) the 
log-divergence seems to be close to reality. 

In order to test the validity of Eq.(24) in the vicinity of crashes, we have 
separated the problems of the divergence itself and the oscillation conver- 
gences on the other hand, in order to extract two values for the rupture point 
tc- (i) tcdiv for the power (or logarithmic) divergence and (ii) tcosc for the 
oscillation convergence. In so doing the long range and short range fluctua- 
tion scales are examined on an equal footing. The final tc is obtained at the 
intersection of two straight lines, by successive iteration fits. The results of 
the fit as well as the correlation fitting factor R have been given in [ |66p7[ |. 

A technical point is in order : The rupture point tcosc is estimated by 
selecting the maxima and the minima of the oscillations through a double 
envelope technique Finally notice that the oscillation basic frequency 
depends on the connectivity of the underlying space • The log-periodic 

behavior also corresponds to a complex fractal dimension p2| , p8[ . 
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